One-Electron Spectral Functions of 
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We calculate the one-electron spectral function of the attractive (negative- (7) Hubbard model. 
We work in the intermediate coupling and low density regime and obtain the self-energy in an 
approximate analytical form. The excitation spectrum is found to consist of three branches. The 
results are obtained in a framework, based on the self-consistent T-matrix approximation, which is 
compatible with the Mermin- Wagner theorem. 
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(A) INTRODUCTION 

The Hubbard model involving electrons on a lattice, 
subject to an attractive interaction when they are on the 
same site, is one of the simplest models for describing su- 
perconductivity. Despite of its simplicity, it has turned 
out to be very challenging for the theoreticians to give 
a simple description of its properties which is valid in 
the various regimes of coupling strength. In the weak 
coupling regime, the link with BCS theory of supercon- 
ductivity has been done by Nozieres and Schmitt-Rink 
0. At sufficiently low T, an instability of the Fermi sea 
towards superconductivity occurs. In three dimensions, 
the transition is essentially mean-field in character. In 
the opposite strong coupling hmit {\U\ — > oo), the elec- 
trons form bound pairs which are immobile since they 
can only move via virtual ionization with an infinite en- 
ergy barrier. However, for large but finite U , those bound 
pairs essentially behave like heavy hard core bosons (with 



an effective mass m* 



-) which are undergoing Bose- 



Einstein condensation at sufficiently low T. On the lat- 
tice, Tc vanishes in the limit \U\ — > oo, while in the 
continuum limit, it remains finite [Q. This difference is 
due to the absence of a pair hopping term when working 
on the lattice. 

In the intermediate coupling regime, the physics will be 
dominated by the interplay between the quasi-particles 
and the bound pairs, which may lead to non-trivial be- 
havior. 

Even though it is still lacking a microscopical deriva- 
tion, the model is interesting in its own right, since it 
allows for studying various routes leading to supercon- 
ductivity. Since the interaction is local, it will be s- 
wave superconductivity, but the generalization to non- 
local interaction can be considered 1^. In the weak cou- 
pling regime, perturbation theory is expected to work, 
and this has been worked out by a series of authors 
0, some of them focusing on 2-d systems 1^. Of spe- 
cial interest is the low-density regime where chances of 
obtaining meaningful results are better, since the ra- 



tio of the scattering length to the average inter-particle 
distance can be used as a small parameter. Unfortu- 
nately those calculations are quickly becoming very in- 
volved since the simplest conserving approximation is 
the self-consistent T-matrix approach [^[0,^. Alter- 
natively, Variational-Monte-Carlo (VMC) calculations, 
based on the Gutzwiller wave-function ||^ and Quantum- 
Monte-Carlo (QMC) simulations have been performed 
JlOf . These methods are providing results which then 
generate a need for a qualitative analytical understand- 
ing. To that aim simpler calculations based on Hubbard- 
Stratonovitch decoupling of the interaction |11|, slave- 



boson mean-field calculations (see for instance Ref. |12| ] 
and references therein), or on the moment calculation of 
the electronic spectral function have been performed . 
Unfortunately the latter does not account for the damp- 
ing of the quasi-particles. 

The aim of this paper is to treat analytically the interme- 
diate coupling regime, which is the most delicate. This 
allows us to give an analytical account of the results ob- 
tained with QMC simulations. We first review the self- 
consistent T-matrix approximation. As pointed out by 
several authors fl^ , the corresponding numerical calcu- 
lations typically yield a superconducting instability at a 
finite T, even in two dimensions. This contradicts the 
Mermin- Wagner theorem. We then propose an alterna- 
tive scheme which complies with this theorem. We then 
proceed to the calculation of the electronic structure. 



(B) THEORETICAL FRAMEWORK 

We study the Hubbard model on the square lattice : 
H = Yl^^ioAaC].'y + U^n,^^ni^i. (1) 

i,j (a) i 

We consider an attractive interaction {U < 0) in the 
intermediate coupling regime {\U\ < W), W being the 
band width. In 2 dimensions, any attractive potential 



1 



has a bound state. In the case \U\ = W, the binding 
energy Eb has been found to be Eb « 0.2W 0, namely 
Eb <C W. We are thus in a situation where bound pairs 
exist and have a strong influence on the physics via the 
splitting of the non-interacting band into 2 sub-bands. In 
this regime, the pairs are extended. They become purely 
local only in the \U\ = 00 limit since for any finite U 
they can move via virtual ionization j^]. We also note 
that the BCS theory successfully describes the weak cou- 
pling regime. However there does not exist any analytical 
theory in the intermediate coupling regime, and most re- 
sults are obtained out of numerical simulations HJl^ . In 
the low-density regime, the self-consistent T-matrix ap- 
proximation is expected to be exact and has been solved 
by a variety of authors Unfortunately, numer- 

ical difficulties prevented those authors from obtaining 
results for arbitrary U. We also note that the numerical 
solutions may lead to unphysical results such as a finite 
critical temperature for Bose condensation of the pairs 
in two dimensions, which is contradicting the Mermin- 
Wagner theorem. We believe (see below) that this is due 
to the use of an inappropriate expression for the particle 
density. That however does not discredit the scheme, and 
we are basing our approach on it. It amounts to solving 



-U 



1 + Ux{q,ii^n) 



(2) 



X{q, ivn) = 



G{p, iuJn)G{q- p, iVn - iujn) (3) 



iujn - + ^ - T,{q, iujn) 



(5) 



Here, w„ are Fermionic, and i^n Bosonic Matsubara fre- 
quencies. This set of equations is valid above Tc, as no 
anomalous Green's function enters. Otherwise one can 
resort to the scheme obtained by Pedersen et at, by func- 
tional derivative techniques This approximation is 
conserving and it diagrammatically corresponds to sum- 
ming up the dressed particle-particle ladder which in- 
cludes the leading order in an expansion in kpa 
Another important quantity is the two-particle Green's 
function which is defined by : 

G(2) {q, lu^) = e-"- (T, [Q(q, r)gt(-q, 0)] ) dr, (6) 
JO 

where T,- is the usual time-ordering operator and the op- 
erator 



(7) 



creates a pair having (center-of-mass) wave vector q. 
G^^^ (q, iVn) is related to the T-matrix by : 



U + T{q, iVn) 
C/2 



(8) 



We calculate G'''^\q,ivn) by inserting the free-electron 
Green's function into expression (||) for xiqJ'^n)- For 
simplicity, we approximate the density of states (DOS), 
p{e), of the tight-binding band resulting from the Hamil- 
tonian (|l]) by the square DOS (i.e. p{e) — ^ for |e| < ^ 
and p{e) = otherwise). 
For small momenta, G^'^^{q, ivn) is given by : 

G^''Hq,lUn) = 



21^(1-1, 



ln$- l + In 



iVn+tlB — \Eb\- 



(9) 



where ^5 2/i + T4^ + \Ebl $ = \^b\ is the 

binding energy of a pair. The binding energy is obtained 
as a solution of : 



1 

U 



x{q = o,t^ = £^6)lu=- 



which yields 



\Eb\ = 2W 



(10) 



(11) 



The form (||) has the correct behavior for Vn going to in- 
finity in the low density regime, i.e. G^'^^ {q, ii'n) 

The spectrum of G'-^-' ((f, w„) presents two features : (i) 
a sharp quasi-particle peak, which can be found by ex- 
panding (|^) with respect to ii^n + ^ (ii) ^ continu- 
ous spectrum which extends over energies above the one 
of the quasi-particle. Correspondingly, the lowest order 
form of the T-matrix, valid for small wave vector and 
frequency is given by : 



Toiqjiyn) 



(12) 



The mass renormalization factor of a pair is given by : 



m 

TO," 







W+\Eb 
W 



2H/2 



<f>ln<f>. 



(13) 



In the intermediate coupling regime, the mass is only 
weakly rcnormalized while in the strong coupling regime, 
there is a strong renormalization of the order . Due to 
the relationship between the two-particle T-matrix and 
the two-particle Green's function, the quantity /is that 
we defined above does represent the chemical potential 



2 



of a pair, which has bosonic character. 

For q's close to the nesting vector Q = (tt, tt), we obtain: 



4m* 



2^l+\u\ 



(14) 



In the vicinity of the zone corner, the renormaHzation of 
the pair-mass is different from the one close to the zone 
center. Even in the intermediate coupling regime, it is 
strongly renormalized to be ^ ~ At g = Q, the 
form (9) of the T-matrix is actually exact, related to the 
fact that the creation operator 



^ P+Q. 



(15) 



of an "77-pair" with center of mass momentum Q, satisfies 
the simple commutation relation ||l^,|l^,|l^ : 



(16) 



Using the above found expressions ([l^ ) and (|lj), we can 
calculate the self-energy. To lowest order, we insert the 
free-electron Green's function in Eqn. (^). The first 
contribution to the self-energy arises from the poles of 
the T-matrix. Due to the statistical factors wc obtain 
(to that order of approximation) that the contribution of 
the 77-resonance is exponentially small, as well as those 
following from the poles of the Green's function. After 
performing analytical continuation, we are left with : 



Si(fc, w) 



uj + tj: - H + Hb + ?'0^ 



(17) 



The quantity will be defined below, in Eqns. <J^^ and 
(^). Ei(fc, w) yields then the Green's function as : 



G{k,uj) = 
2 I xr. 



1 f _2jt^-j^)j^ 



1 



+ 5MS — ^xj: + iQ+ 
1 



(18) 



where x-j: — 
diately note 
momentum k 



imme- 
espect to 



J {2 {tj: - fi) + fiB^ + AU^fid. We i 
the two limiting behaviors, with respi 

Xi:Ki A + 2t-fk^, (19) 



with 7 — and 



for small momenta; respectively. 



(20) 



2^fc « 2 iu -n)+^lB 



2U^nd 



2 (tz ~ fi) + (Ib' 



(21) 



for large momenta. 

At this stage of the calculation, the Green's function has 
a two-pole structure. The lower excitation branch cor- 
responds to quasi-bound fcrmions (hereafter denoted as 
"bosonic" band), while the upper branch describes the 
unpaired fermions (fcrmionic band). At small momenta, 
we obtain: 



G{k,uj) 



2A w + i (/is - A) - 7t/c2 + io+ 

A-|gfe| 1 

2A cj-fi(^B+A)+7ifc2 + iO+ 



(22) 



with the spectral weight mainly located in the unpaired 
fermion band (first contribution in Eqn. (p^)). At large 
momenta, the Green's function results into : 



1 



G{k,uj) 



+ 



(2(*£-A^)+A'fl) 

oj- {tj:~ 11)+ iO+ 

(2(t£-^')+A'B)' 

u> + tj: - ^ + ^B + 



(23) 



where the weight of the paired fermion band is even 
smaller than for small momenta. The form of the Green's 
function Eqn. (|lj) differs from the one of Ref. 11] 
because the chemical potential fi is located below the 
fermionic band in our problem. 

We note that there are two equivalent expressions for the 
particle density operator h: 

fli = ^ (?li,o-(l - fii -cr) + fli^^fli^^a) (24) 
a 

On one hand, we can use the left hand side to express 
the particle density n as ni, where the subscript 1 indi- 
cates that the density is calculated out of the one-particle 
Green's function : 



(25) 



iuJn ,k 



Alternatively, we may use the r.h.s. of Eqn. ( p^ ) by 
separating explicitly the contributions from the unpaired 
fermions (first term) and the doubly occupied sites (sec- 
ond term). The total density rid of the latter is given 

by 



E r -Irn 

4- ./-^ 2^ 



{G(2)(q,co + iO+)}iVB(a'), (26) 
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Nb{i^) being the Bose- Einstein distribution function. 
Owing to the latter, at low temperatures, only 
the low-energy part of the two-particle spectrum, 
Im {G^^^((7', w -I- «0+)}, will contribute to n^. According 
to the discussion following Eqn. (pl|), this low-energy 
part has two contributions : the sharp quasi-particle ex- 
citation, given by (p^), which dominates for small mo- 
menta, and the low- frequency tail of the continuous spec- 
trum of G^^^ (g, ii^n)- Since it does not have a sharp struc- 
ture, it would only contribute a featureless background 
to the one-electron spectral function; thus, we have ne- 
glected it for the calculation of the self-energy 'Si{k,uj). 
Neglecting the continuum also in calculating n^, we find, 
according to Eqns. (|), (|l|) and (||) : 



nd 



Am 



-Z 



(27) 



Thus, in our approximation, the number of doubly occu- 
pied sites Ud is given by the "number of bosons" , hb 



TIB = ^ iVs 



weighted by a factor 



4rn 



fJ-B 



\2W ) ^ ' ' 



(28) 



(29) 



which is the residue of the two-particle Green's function 
at the bottom of the two-particle band. These bosons, 
having energy Eb — and chemical potential /i^, 

represent pairs of electrons being (virtually) bound by the 
on-site attraction. These results correspond to the obser- 
vation of other authors (see, for example Ref. [|j) that, 
for sufficiently strong attraction, the two-particle Green's 
function, (respectively the T-matrix) can be interpreted 
as a "bosonic Green's function". However, there is a 
weight factor between the two. For large \U\, this weight 
factor R goes to one : in the strong coupling limit, all 
the double occupancy is due to coherently propagating 
quasi-bound pairs. In the intermediate coupling regime, 
which we want to consider {\U\ < W), the weight factor 
R is roughly 0.5 : only about one half of the two-particle 
spectrum in Eqn. (^ is resulting from coherent exci- 
tations. In the weak coupling regime, R vanishes, thus 
rendering our approximation invalid in this limit. 
Since, in our approach, there is gap in the one-electron 
spectrum, separating the "bosonic" band from the 
fcrmionic one, there is no difficulty to obtain the number 
of unpaired electrons np as : 



(30) 



where /F(eg^) is the usual Fermi distribution function. 
The dispersion ^ of the unpaired fermions and the 
spectral weight entering Eqn. ( |30| ) are : 




The total particle density results into 



(31) 



(32) 



For non- interacting particles, we may equally well use 
both r.h.s. and l.h.s of Eqn. ( p^ to calculate the density, 
since Wick's theorem applies. However, the identity (|24| ) 
may be violated in an approximate treatment like pertur- 
bation theory. This is the reason why the self-consistent 
T-matrix calculation breaks down at low T in two di- 
mensions. By calculating the density by means of Eqn. 
(^), there is nothing preventing T{q — Q,uj — Q) from 
diverging at finite T , signaling a phase transition, while 
using Eqn. ( ^2|) would definitively keep r(0, 0) finite for 
any finite temperature. By making use of Eqn. (^2|), we 
make sure that Bose condensation can only take place at 
T = 0, in agreement with the Mermin- Wagner theorem. 
Indeed, according to the expressions ( |27| ) and ( p8| ) for Ud 
and ns, respectively, the bosonic chemical potential /i^, 
for a given Ud (or ub) and for d = 2, is different from 
zero at any finite temperature, which inhibits Bose con- 
densation, except for T = 0. In principle, for d = 2 one 
should see a Kosterlitz-Thouless (KT) transition. How- 
ever, since our approximation does not treat the bosonic 
phase fluctuations in an adequate way, we cannot expect 
to see the KT-scenario. On the other hand, Tc may well 
be finite for a 3-d system. Actually our procedure is sim- 
ilar in spirit to the two-particle self-consistent approach 
to the repulsive Hubbard model by Vilk et al. p8| , p^ . 
Now for a two-dimensional system, we can explicitly eval- 
uate the number of pairs. We obtain : 



rid 



np 



2R 

pwz 



In 



1 



1 



\Eh 



In 



62 



1 



(33) 



(34) 



62 



where A and 7 have been given above. 
We are now ready to (numerically) solve Eqns. (|3^), 
( ^ ) and (IIJ) for the chemical potential, as a function 
of temperature. Hereafter, the density n is fixed to be 
n = 0.1. The resulting ij,b{T) vanishes exponentially at 
T = 0. From Eqns. (|) and (^, we obtain that the 
range of the two-particle Green's function, ^, is given by 



t 



2|Mi3|' 



(35) 



where a is the lattice spacing. It is displayed on fig. 1, 
for three different values of the coupling strength U. We 
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note that it shows a strong dependence on U . At very 
low temperatures, it is given by : 



Bt llWZr, 

— e 
2 



(36) 



Oppositely, in the intermediate temperature range, ^ be- 
comes independent of U. We also find that ^ decreases 
with increasing \ U\. In the strong coupling limit the ratio 
^ tends to zero and thus the exponential divergence of 
^ in y is suppressed (there remains only the power-law 
dependence ^ ~ '^)- 

We may now define a coherence temperature Tcoh as 
the temperature at which the range of the two-particle 
Green's function exceeds 10 lattice spacings. We obtain 
Tcoh « 0.16t for U = -U and T^oh « O.li for U = -6t, 
which may be compared to Tc as obtained from the nu- 
merical simulations . We see that they compare favor- 
ably and that, moreover, Tcoh decreases with increasing 
U. Finally, we note that ^ becomes of the order of the 
lattice spacing at temperatures well below \Ei,\/2. We 
now turn to the temperature dependence of the number 
of pairs. It is displayed on fig. 2, as a function of T, for 
three different values of U. At low temperatures, it is 
independent of T. Since the binding energy of the pair 
(and the gap) decreases as U gets smaller, Ud begins to 
decrease at a lower T for J7 = — 6^ than for U ~ —lOt, 
for example. In all cases decreases by a factor 2 at 
T « A/2. 

In order to reach a better self-consistency, we now cal- 
culate how the quasi-bound states affect the two-particle 
propagator. Introducing Gi(g, w) as the lower branch of 
G{q,uj) given by Eqn. (|l^), we can calculate the first 
correction to the two-particle propagator as xi, with : 

+ ^ Go((f- p, iVn - iuJn)Gi{p, iujn). (37) 

Carrying out the summation over Matsubara frequencies, 
we obtain that the statistical factors are exponentially 
small and thus there is no correction to x to that or- 
der. Thus, the T-matrix is still given by Eqn. (p^. We 
can now proceed to the calculation of the second-order 
correction to the self-energy. It is given by : 

S]2(g, iujn) = ^ To(p, Wn)Gi{p + q, iVn + ii^n), (38) 
yielding, after performing analytical continuation, 

^2{q,oj) = 



Xq - 2 (t^ - /i) - flB 



uj - ^{x^- 11b) + iO+ 



W 



(39) 



where ujc is a frequency cut-off needed by the assumption 
that the 77-resonance is sharp for ^'^ ^} < ujc ^ ^ , with 



TO* = TO-Lp. This corresponds to |Q — ifl < as it is 
borne out by numerical calculations of the two-particle 
Green's function G^"^^ covering the whole Brillouin zone 
pT| , p7t . The precise value of the cut-off has little infiuence 
in the following (numerical) results. We note that the 
first contribution in Eqn. ( |39| ) is following from the long- 
wavelength behavior of the T-matrix, while the second 
is due to the ry-resonance. Even though both Si and 
the first contribution to S2 originate from the pole of 
To((f, Wn), given in Eqn. (|l^), they are found to have 
opposite dispersions. The total self-energy results as : 



I](g,tj) = Si((f,cj) + ^2{q,i^)- 



(40) 



At this point of the calculation, we can recalculate the 
particle density ni by evaluating explicitly Eqn. (^5|), us- 
ing the Green's function resulting from Eqn. (|4^) . Since 
there is no a priori reason that this would yield a result 
comparable to what is following from Eqn. (^), this is 
a consistency check of the framework we are using. The 
result is displayed on fig. 3. It is obvious that this is in 
very good agreement with the expected value n = 0.1 for 
all values of U. This means that our calculation is con- 
sistent. It also implies that, on one hand a small change 
in the chemical potential corresponds to a small change 
in the particle density, but on the other hand, it may 
induce a big change in the pair density. 



(C) EXCITATION SPECTRUM 

We can now summarize our findings by plotting the 
various pieces of the spectral function. This is done on 
fig. 4a, for U = ~6t {Eb = -t) and on fig. 4b for 
U = —lot {Eb — — 4.2<) for momenta along the diagonal 
of the Brillouin zone. In both cases, the spectrum con- 
sists of four branches. The spectral weights of the various 
branches are displayed on fig. 5a for U = —6t and fig. 
5b for U = -lot. 

These results can be commented as follows : 

(i) There is a branch at negative energies following from 
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the quasi-bound states. It is centered around — 
with a width W - < w. The de- 

pendence in U is weak, and the dispersions for \U\ — 6t 
and \U\ = lot are essentially the same, up to a small 
shift; moreover, we note that they are opposite to the free 
fermion dispersion. The weight decreases with increasing 
momentum and gets mostly negligible for ^ 2. It de- 
creases with increasing \U\, at small q 's, but the total 
weight remains constant. 

(ii) The fermionic band that was coming out of the first 
approximation (Eqn. (|l^)) is now resulting as a super- 
position of two branches which would merge into a single 
one, were damping taken into account. This superpo- 
sition is done out of the two branches which are most 
free-electron-like. This is somewhat arbitrary, especially 
in the domain where the hybridization with the other 
branch, which lies at positive energies too (see below), 
is strong. Consequently, we do not show the spectral 
weights in this range. The width of the fermionic band 
is slightly larger than the original bandwidth, and shows 
little dependence on. U. It lies at slightly larger energies 
for \U\ = lot than for \U\ — 6t, but again the dependence 
in U is weak. Its dispersion is parallel to the free fermion 
one. 

(iii) There is a third branch resulting from the rj- reso- 
nance; its dependence in U is weak and both curves (for 
\U\ = 6t and \U\ — lOt) are parallel. It comes down 
in energy with increasing \U\. The contribution of the 
ry-resonance band is very small, barely accounting for 5 
percent of the spectral weight at its maximum, which is 
located at q = Q. 

It is interesting to compare our results with Quantum 
Monte Carlo (QMC) simulations [|o|j20|. The calcula- 
tions of Ref. have been performed for n — 0.4 and 
for \U\ = -At, -8t, -12t. This is hkely to be outside the 
realm of densities for which the T-matrix approximation 
is really valid. Nevertheless, the following features coin- 
cide with our findings : there are two distinct excitation 
branches which get more and more clearly separated by 
a gap when \U\ increases. The width of the two bands 
is smaller than in our calculation, in particular the lower 
(bosonic) band has very little dispersion. However, the 
weights (the fermionic band has large weight for large 
wave vectors, whereas the weight of the negative energy 
bosonic band is concentrated near the zone center) is sim- 
ilar to our figures 4 and 5. 

Ref. presents new results for the one-electron spectral 
function for n = 0.1, the same value that we have used, 
and \U\ = 8t. Here again, there is a fermionic branch at 
positive energies E (positive with respect to the chemical 
potential), separated by a "pseudo-gap" from the excita- 
tions at ii^ < 0. Weight and width of the fermionic band 
is similar to our result. The excitations at i? < have 
most weight near the zone center, as we find it. However, 
their structure seems to be more complex : with increas- 



ing wave number they split into two " sub-branches" . The 
lower part has downward dispersion and seems to cor- 
respond to our bosonic branch. The other sub-branch 
produces weight near the chemical {E ~ 0) potential for 
large wave numbers. This might correspond to our 77- 
resonance, provided that its energy near the zone corner, 
lying at a positive energy in our calculation, would in 
reality be lower. In this respect, the spectral functions 
for n = 0.4 and \U\ — ~6t, also shown in Ref. |2^, are 
particularly interesting : ioi E > 0, besides the strong 
fermionic branch, there is a second branch of excitations 
the dispersions of which is very similar to the ry-peak in 
our figures 4 and 5. 

Thus, the main features of our spectral functions seem to 
be present in the Monte Carlo results, although the lat- 
ter show a more complex structure. A detailed and more 
quantitative comparison with QMC should also take into 
account the fact the "Maximum Entropy method" used 
there in order to extract spectral functions from data 
obtained as functions of the (imaginary) Matsubara fre- 
quencies does not easily allow for an unambiguous iden- 
tification of excitations with small weight. 
Finally, we note that we do not obtain any spectral 
weight at zero frequency, signaling the presence of a cor- 
relation induced gap. We also checked that including 
particle-hole excitations in the calculation does not af- 
fect this conclusion. The appearance of a true gap in our 
calculation is, at least partly, due to the fact that our 
spectral lines have no width (except that the " doubling" 
of the fermionic branch gives a hint to a broadening of the 
latter). Spectral functions with finite line width would 
be obtained either by doing the g-sums in the expression 
(^) for the self-energy more precisely, and/or by evaluat- 
ing G and E by solving Eqns. (|^) to (^ self-consistently. 

In summary, we have determined the excitation spec- 
trum of the attractive Hubbard model at intermediate 
coupling out of a simple analytical calculation. We first 
pointed out an intrinsic problem of perturbation theory 
relative to the implementation of the Mermin- Wagner 
theorem. We made use of an alternative expression for 
the density to obtain a qualitatively correct theory which 
does not break down at low T in two dimensions, in con- 
trast to previous self-consistent calculations. We obtain 
an analytical expression for the Green's function which 
reproduces the qualitative features of the QMC simula- 
tions in the low density regime. 
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FIG. 4. a) Spectrum of the Green's function at n = 0.1, 
T/t — 0.1 and U = — 6t It consists of the bosonic band 
(dashed line) , the fermionic bands (dotted and dashed-dotted 
lines) and the r;-resonance band (full line), b) The same for 

U = -lot. 



FIG. 5. a) Spectral weight of the bosonic band (dashed 
line), the fermionic bands (dotted line) and the T^-resonance 
band (full line), for n = 0.1, T/t = 0.1 and U = ~6t. b) The 
same for U — — lOt. 



(1996) and |cond-mat/970218S . 

J.M. Singer and T. Schneider, private communication. 



FIG. 1. Range of the two-particle Green's function, as a 
function of temperature T, at density n = 0.1, for three dif- 
ferent interaction strengths, U = — 4t (short-dashed line), 
U = -6t (full line), U = -8t (dotted line) and U = -lOt 
(long-dashed line). 



FIG. 2. Density of quasi-bound pairs as a function of tem- 
perature T, at density n = 0.1, for U = —6t (full line), 
U = ~-8t (dotted line) and U = - Wt (dashed line). 



FIG. 3. Particle density ni which is obtained by using Eqn. 
(25); the chemical potential is calculated numerically out of 
Eqn. (32). 
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